OU-HET- 740/2012; MISC-2012-04 



en 

1) 



Abstract 



Yukawaon Model with U(3)xS3 Family Symmetries 

Yoshio Koide a and Hiroyuki Nishiura 6 

a Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 
E-mail address: koide@het.phys.sci.osaka-u.ac.jp 
b Faculty of Information Science and Technology, Osaka Institute of Technology, Hirakata, 

Osaka 573-0196, Japan 

(N 1 

E-mail address: nishiura@is.oit.ac.jp 

(N 

A new yukawaon model is investigated under a family symmetry U(3) x S3. In this model, 
*-£h I all vacuum expectation values (VEVs) of the yukawaons, (Yf), arc described in terms of a 

fundamental VEV matrix (<&o) as m the previous yukawaon model, but the assignments of 
quantum number for fields are different from the previous ones: the fundamental yukawaon 
$0 is assigned to (3,3) of U(3)xU(3), which is broken into (3, 1 + 2) of U(3)xS 3 , although 
quarks and leptons are still assigned to triplets of U(3) and yukawaons Yf are assigned to 
6* of U(3). Then, VEV relations among Yukawaons become more concise considerably than 
" the previous yukawaon models. By adjusting parameters, we can fit not only quark mixing 

parameters but also lepton mixing parameters together with their mass ratios. 
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1 Introduction 

It is interesting to consider that the observed quark and lepton mass spectra and mixings 
can be understood from a "family symmetry". In the standard model (SM) of quarks and 
leptons, their mass spectra and mixings originate in the structures of the Yukawa coupling 
constants. If we intend to understand the observed mass spectra and mixings by adopting a 
non-Abelian gauge symmetry as a "family symmetry", the symmetry will be explicitly broken 
at the beginning because the Yukawa coupling constants have family indices. Therefore, most 
family symmetry models are based on a discrete symmetry. However, there is an easy way to 
consider non-Abelian gauge symmetry as the family symmetry: the Yukawa coupling constants 
are effective ones (/ = u, d, e, • • • ) and those are given by vacuum expectation values 

(VEVs) of new scalars Yf. 

(Yf f h = f <(r,)y>. (1.1) 

The fields Yf are called as "yukawaon" pQ. The yukawaon model is a kind of "flavon model" 
[2], but, differently from the conventional navon model, the effective Yukawa coupling constant 
Yj^ is essentially given by a 3 x 3 VEV matrix of a single scalar Yf. Here we do not consider 
that Yy J is given by a linear combination of scalars Ya, Yb, • • • , which belong to different 

representations of a family symmetry G, as Yf = c\(Ya) + C2(Yb) ~\ ■ 

In the previous yukawaon model, the most characteristic feature is that all VEVs (Yf) are 
described in terms of only one fundamental VEV matrix (3> e )- For example, mass matrices for 



quarks and charged leptons are given as follows [3]: 



M e oc(Y e )oc($ e )($ e >, 

M d oc (Y d ) oc {$ e )(l + a d X){<S> e ), 

mU 2 oc (* u ) oc ($e)(l + a u X)($ e ) 



(1.2) 
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For neutrino sector, we consider a seesaw type of neutrino mass generation M v = m^M^ mj^, 
where the Dirac neutrino mass matrix mo and the right-handed neutrino Major ana mass matrix 
Mr are given, respectively, by tod oc M e and 



where {P u ) = diag(+l, — 1, +1) in the diagonal basis of the up-quark mass matrix M u . The 
model [3] has roughly well described masses and mixings of the quarks and the leptons with 
a few parameters. However, the fitting of the quark mixing has been somewhat unsatisfactory 
compared with those of the neutrino mixing. (For more precise parameter fitting, see Ref.[3].) 

In the present paper, we propose a new yukawaon model where all yukawaon VEVs are 
described in terms of a new fundamental VEV matrix (<3?o) similar to the previous model. 
However, in contrast to the previous model (1.2), the charged lepton mass matrix M e is given 
by the same structure as those in the down- and up-quark sectors: 



[For more details, see Eqs.(4.9) - (4.14).] In the past yukawaon models, it has been assumed 
that the factor (1 + a q X) (q = u,d) in Eq.(1.2) originates in VEV structures of additional 
fields. However, in the present paper, we assign the fundamental yukawaon $o to (1 + 2) of a 
permutation symmetry S3 [5], and we consider that the factor (1 + a q X) is due to a coefficient 
of (1 + 2) x (1 + 2) — > 1 under S3. Therefore, it can naturally be understood that all the factors 
which are sandwiched by ($0) are given by (1 + dfX) as shown in Eq.(1.5), instead of Eq.(1.2). 
The details will be discussed in the next section. 

In Sec. 3, we investigate possible superpotential forms for yukawaons, and in Sec. 4, we 
summarize VEV relations in the present model. In Sec. 5, we give parameter fittings for the 
observed lepton mixing (PMNS mixing [6]) and the observed quark mixing (CKM mixing [7]) 



M R oc (Y r ) oc (Y e )(P u )($ u ) + {<S> u ){P u )(Y e } + ■■■ 



(1.4) 



M e oc (Y e ) oc <$ )(l + a e X)<$o>, 
M d oc (Y d ) oc <<D >(l + a d V)(<D >, 
Ml' 2 oc {%) oc <$ )(1 + a u X)(^ ) 



(1.5) 



together with their mass ratios. In the present model, we have 8 adjustable parameters a e , a u e iau , 
a d, and (01, ^2) for 16 observables (2 up-quark mass ratios, 2 down-quark mass ratios, 2 

neutrino mass ratios, 4 CKM mixing matrix parameters, 4+2 PMNS matrix parameters). The 
final section 6 is devoted to summary and discussions. 

2 Fundamental Yukawaon $ and a permutation symmetry S3 

First, we give a brief review of the factor (1 + dfX) in S3 symmetry. When we denote a 
doublet (ipirjipn) and a singlet ip a in a permutation symmetry S3 [5] as 




75(^1-^2) 
^(^1+^2-2^3) 



(2.1) 



tpa = +^2+^3), 



(2.2) 



the field ijj = (V'i>V'2>V'3) i s represented as 
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Then, a bilinear form ipi/j is invariant under the S3 symmetry only when tpa^ab^b is given by the 
form 

IpaZablpb = ^a(l + a f X) ab 1p b , (2.4) 

where aj is a free parameter, and 1 and X are defined by Eq.(1.3). The appearance of a free 
parameter a/ is due to a reason that there are two singlets which are composed of tp a , i.e. ipatpa 
and (tp 7V ,tp rl ) T (i/} 7r ,ip rj ). 

As seen from Eq.(2.4), we can assume that the fundamental yukawaon $0 is transformed 
as (3,2 + 1) of U(3)xS 3 , i.e. 

^(l + a/lf^. (2.5) 

However, since the bilinear form (2.5) cannot have quantum numbers which distinguish the 
sectors f = u,d,e {at are merely free parameters), we modify the expression (2.5) into 



(2.6) 



where Sf are fields and indices a,/3 are of another U(3) symmetry (we denote it as U(3)'). We 
assume that U(3)' is broken into S3 at an energy scale of A' by non-vanishing VEV {Sf) whose 
forms are given by 

(S f )=v Sf (l + a f X). (2.7) 

Here, we assume A' S> A (U(3) family symmetry is broken at A). Note that the indices a,f3 in 
the expression (2.6) are of U(3)', while the indices a, b in the expression 



(*o)£<S?>(* )y 



(2i 



are of S3. It is worthwhile to notice that in the conventional S3 family model quarks and leptons 
are assigned to (singlet+ doublet) 's of S3, while in the present model quarks and leptons are 
assigned to triplets of the U(3) family symmetry, and only the fundamental yukawaon $0 is 
assigned to (3, 1 + 2) of U(3)xS 3 . 

3 Superpotential 

In the present model, would-be Yukawa interactions are given as follows: 



W Y 



y, 



^l 5 H d + J^ff, + \rvIY 1 £v] + fdfflqjHt + ^Y^ qj H u , (3.1) 



Vv 



ydic v ij, 



where I = (vl-,&l) and q = (ul,cIl) are SU(2)l doublets. Under this definition of Yf, the 
CKM mixing matrix and the PMNS mixing matrix are given by Vqkm = utll d and Upmns = 
UeU u , respectively, where Uf are defined by UjMjMfUf = (Df are diagonal). In order to 
distinguish each yukawaon from others, we assume that Yf have different R charges from each 
other together with R charge conservation. (Of course, the R charge conservation is broken at 
the energy scale A'.) Here, we have assumed that the R charge of the Dirac neutrino yukawaon 
Y v is identical with that of the charged lepton yukawaon Y e , so that Y v is replaced with Y e for 
an efficient use of fields. 

We assume the following superpotential for yukawaons: 
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(3.3) 



W u = - 



KEikYu l E[j + \' u Y£E Kl Yij 



A'' 



n' u Y% + f(<S>o)lSf(S> ) bj 
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W R = j {\ R E lk Y^E l3 + \' R [y&E"^ + Yf k E kl Y k ) + <Tr[Y u E]Y^ + Ti[EY e ]Y^] } eg. 

(3.5) 

(In this paper, we use Y u instead of <E> U unlike the previous papers, since we treat the up- and 
down-quark sectors in the same way Y u <H> Y d .) Here we have assumed that the R charges 
satisfy the following relation 

R(Y e ) - R(Y d ) = R(Y e ) - R(Y d ) = R{S e ) - R(S d ). (3.6) 

The VEVs of the introduced fields E, E and P are described by the following superpotential 



which leads to 



W E ,p = ^-Tr[EEPP] + ^-Ti[EE]Tr[PP], 



ocl, (P)(P)ocl. 



(3.7) 
(3.8) 
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Table 1: Quantum numbers of the fields 
We take specific solutions of Eq.(3.8): 

-(E) = ±(E) = 1, 

VE VE 

— (P) = ^(P) t = diag(e-^,e-^,l), (3.10) 
Vp v p 

as the explicit forms of (E), (E) and (P). Here, we have assigned R charges for those fields as 

R(E) + R(E) + R(P) + R(P) = 2. (3.11) 

Since we consider R(E) + R(E) / R(P) + R(P), terms Tr[EEEE], Tr[PPPP], and so on are 
forbidden by the R charge conservation. However, the absence of the term Tr[EP]Ti[PE] in 
Eq.(3.7) must be assumed ad hoc. 

In the superpotential Wd, Eq.(3.3), we have taken 

R(E) = R(Y d ), (3.12) 

and we have added the m^Eij term to Yf-. The relation (3.12) has been ad hoc assumed in order 
to adjust the quark mass ratio m^i/m^- Such the assignment is not contradict with whole R 
charge assignments of the fields. However, once the assignment (3.12) is done, then, we cannot 
add E term to Y e and/or Y u because of the R charge conservation. 

We list whole fields in the present model in Table 1. As seen in Table 1, the sum of the 
anomaly coefficients A is ^ A = —9. In order to be ^ A = 0, we need further three fields, T^ a , 
Tf a and Tf a . For the time being, we do not specify the roles of these fields in this model. 

4 VEV relations 

We assume that our vacuum always takes (9 a) = (A = e,u,d, •••)■ Therefore, VEV 
relations are obtained from the SUSY vacuum conditions dW/d&A = 0. Since relations from 
another SUSY vacuum conditions dW/dY e = and so on always include one of (9a), such the 
relations do not have meaning except for Eq.(3.7). Since we assume that the SUSY breaking 



(3.9) 



is caused by gauge mediation (except for family gauge symmetries), we consider that our VEV 
relations in the yukawaon sector are satisfied until a low energy scale. 

From the superpotential terms (3.2)-(3.5), we obtain the following VEV relations among 
the yukawaons: 

(Y: j ) = -^ K (P ik )(Y k e l )(P lj ), (4.1) 
<*3> = -^<(*o)£><S?><(*o) W >, (4.2) 
( Y d) = -^( pik )( Y "l)( plj )^ ( 4 - 3 ) 
<*8> = " A(($o)l)(St)((<S>o) bj ) + (4.4) 

w> = -^((^- l )r fe ™(^ m )«„)((^- 1 )r, (4.5) 

= - ^<(*o)£><S?><(* )«>, (4.6) 



A', 



( 1 &[<y tt >(S>](y^> + Tv[(E)(Y^](Y^))] ((E- 1 ))^, (4.7) 

where (S/), (E), (E) and (P) are given by Eqs.(2.7), (3.9) and (3.10), respectively. We have 
assumed that the VEV matrix (<3?o) is diagonal in the basis in which the VEV matrix (Sf) take 
the form (1 + ajX), i.e. 

(*o) = diag(ui,U2,u 3 ). (4.8) 

In the present model, common coefficients are not important. Therefore, when we omit 
those coefficients, quark and lepton mass matrices are given as follows: 



M e = (P)Me(P), 


(4.9) 




(4.10) 


M d = (P)((tf)(S d )(*o)+m° d l) (P), 


(4.11) 




(4.12) 


M v = M e M R l M^ 


(4.13) 


M e + M e M^ 2 + +iu (TV[MV 2 ]M e + Tr[M e ]M^ . 


(4.14) 



For numerical calculations in the next section, we will use dimensionless expressions (<3?o) = 
diag(xi,x 2 ,x 3 ) with x\ + x\ + x\ = 1, (P) = diag(e i<?il , e i<fe , 1) and (S f ) = l + a f X in Eqs.(4.9) 
- (4.14). (Therefore, in Eq.(4.11) is also a dimensionless parameter.) 

5 Parameter fitting 

For simplicity, we assume that the parameter in the down-quark sector is real as well 
as a e in the charged lepton sector, so that only a u is complex. For the convenience, hereafter, 
we denote a u anew as a u e lolu (a u and a u are real). Then, in this model, we have 8 adjustable 
parameters a e , a u e lotu , o^, £ v , m° d and (fa, fa) for 16 observables [2 up-quark mass ratios, 2 
down-quark mass ratios, 2 neutrino mass ratios, 4 CKM mixing matrix parameters, 4+2 PMNS 
matrix parameters]. Since the observed charged lepton mass ratios are used as input values, we 
do not count them as adjustable parameters. The parameter m® is used only in order to fit the 
down-quark mass ratio md/m s . (In other words, other observables are insensitive to the value of 
m^.) Therefore, if we do not count this parameter m^, we have 7 parameters for 15 observables. 

5.1 PMNS mixing 

In order to predict the neutrino mixing (PMNS mixing) parameters, we have 6 parameters 
[a e , £ u , a u e lolu and (fa, fa)]. At present, we know only 5 observed quantities [2 up-quark mass 
ratios, 1 neutrino mass ratio R v = Am^/Aml, 3 neutrino mixing angles, sin 2 29 a t m , tan 2 9 so i ar 
and sin 2 2^13] among 10 observable quantities [2 up-quark mass ratios, 2 neutrino mass ratios, 
4+2 PMNS mixing parameters]. Therefore, in principle, we cannot determine the parameter 
values from the observed quantities. Instead, we use the observables which are described by 
parameters as few as possible. To start with, let us use up-quark mass ratios rf 2 = \/m u /m c 
and r 23 = \J m c /mt which are described by 3 parameters a e and a u e lolu . The observed values of 
r\ 2 and r 23 are as follows [8]: 

r u l2 = J™» = 0.045j™8, r« 3 = .l^= 0.06 ± 0.005, (5.1) 
V m c V m t 

at [i = mz- We illustrate a relation (a u , a u ) and a e which satisfy r" 2 = 0.045 and r 23 = 0.060 in 
Fig. 1. As seen in Fig. 1, the these parameters are bounded in a u = —(1.3 — 1.8), a u = 0° — 10° 
and a e = — 10. (Since the observed center values have large errors, these bounds should not 
be taken rigidly.) 

Next, we focus on the neutrino mixing sin 2 29 a t m . We find that it is insensitive to a u 
and £„, although it is described by parameters a e , a u e !Q ", £ u and (<pi,4>2) m °ur model. In 
Fig. 2, we show curves in the (a u ,a e ) plane which satisfy sin 2 29 at m = 0.90, 0.92, 0.94 and 0.96, 
respectively. Here, for convenience, we have fixed (c/>i, <fo) as (<fii, $2) = (180°, 180°). (For a case 
of (</>i,</>2) = (0°,0°), we could not find reasonable parameter solutions at all.) Here, we have 
used = 0.0025 and a u = 8° and 10°, and (fa, fa) = (180°, 180°) tentatively. As seen in Fig. 2, 
two curves (a u = 8° and 10°) are almost degenerated. For reference, we also show a curve which 
satisfies reasonable up-quark mass ratios. Since we cannot take a u > — 1.3 as we have shown in 
Fig. 1, we cannot obtain a value more than sin 2 26 a t m = 0.97. If we want to obtain a value of 
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Figure 1: Parameter values (a u ,a u ) and a e which satisfy the observed values of up-quark mass 
ratios \Jm u /m c = 0.045 and y 7 m c /m t = 0.060 




Figure 2: Parameter values (a u , a e ) which satisfy sin 2 29 atm = 0.90, 0.92, 0.94, 0.96, respectively. 
(In the figure, "atm" means sin 2 26 a t m ) For reference, the curve which satisfies up-quark mass 
ratios \Jra u jm c = 0.045 and \J m c /mt = 0.060 are also shown. 

sin 2 29 atm as large as possible, we must take a parameter set 

(a e ,a u ,a u )~ (8, -1.35, ±6°). (5.2) 

For this parameter set (5.2), we can obtain tan 2 9 so i ar ~ 0.5 and R ~ 0.04. More accurate 
predicted values will be given after we discuss the CKM mixing values in the next subsection. 

5.2 CKM mixing 

The 4 observables in the CKM mixing matrix depend on 3 parameters and (0i,02) in 
addition to the parameter values (5.2). Although the parameter can be determined by two 
down-quark mass rations [8] 

rf 2 = ^ = 0.053^ rf 3 = ^ = 0.019^, ( 5 - 3 ) 
we cannot give reasonable fitting for the CKM mixing and the down-quark mass ratios simul- 
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Figure 3: Contour Plots in the ((pi, fo) parameter plane, which are shown by using experimental 
constraints on \V US \ = 0.2252 ± 0.0009, \V cb \ = 0.0406 ± 0.0013, \V ub \ = 0.00389 ± 0.00044, and 
\Vtd\ = 0.0084 ± 0.0006 and by taking some values for atm = sin 2 26 atm and = tan 2 9 so i ar . (a): 
sin 2 26> aim =0.98, (b): sin 2 26> aim =0.96, (c): sin 2 2^ atm =0.94, and (d): sin 2 26» atm =0.92, which are 
shown by dashed lines, and solar= tan 2 9 so i ar which is shown by dotted or dot-dashed line, (e): 
tan 2 Osoiar = 0.495 (dotted line), and (f): tan 2 9 so i ar = 0.515 (dot-dashed line). Here we take 
the values for the parameter set of [a e , {a u ,a u ),S, v ,ad,m d ] given in Eq. (5.10). 

taneously. Therefore, we fit only the value of r| 3 by ad, and the value rf 2 is adjusted by the 
parameter m d which does not affect the CKM mixing parameters. 

Although we have roughly obtained the parameter values of (o e , a u , a u ) in Eq.(5.2) from r 23 , 
r" 2 , and sin 2 2Q atm , the predicted values for CKM mixing are also dependent on the parameter 
values (a e , a u ,a u ). We consider that the observed values for the quark mass ratios are still contro- 
versial. Therefore, we search a set of reasonable parameter values for [a e , (a u , a u ),£ v , cid, (</>i, 
that can give the observed PMNS and CKM mixings within one sigma and quark mass ratios 
within two sigma, keeping the result (5.2) in mind. Namely, After searching a rough value of the 
parameter which can give reasonable CKM mixings and the down-quark mass ratio r^, we do 
a fine tuning of parameter values for [a e , (a u , a u ), a^, (4>i, &)] with help of figures demonstrated 
in Fig. 3. 

As a result, we obtain the following predictions: 

sin 2 20 at m = 0.965, tan 2 B so i ar = 0.522, sin 2 26 13 = 0.027, (5.4) 
S e cp = -177° (J e = -9.3 x 10~ 4 ), (5.5) 

R u = s ° lar = 0.044, (5.6) 

Am atm 

\V US \ = 0.2240, \Vcb\ = 0.0404, \V ub \ = 0.00409, \V td \ = 0.00823, (5.7) 
<*9 p = 65.3° (J g = 3.3 x 10~ 5 ), (5.8) 



r\ 2 = 0.0422, r& = 0.0659, rf 2 = 0.0530, r& = 0.0360, (5.9) 
under the parameter values 

a e = 8.5, a u = —1.32, a u = —6.5°, = 17, m[J = 0.0113, 

£, = 0.0019, (0i, <fo) = (174.5°, 195.9°). (5.10) 

Although the predicted value tan 2 so iar = 0.522 is somewhat large compared with the observed 
value [9] tan 2 so j or = 0.468^qq29> it is consistent with the KamLAND data [10] tan 2 6 so [ ar = 
0.56lQQ7(stat)^QQg(syst). The predicted value sin 2 26^ = 0.027 is small compared with the T2K 
data [11] 0.03 < sin 2 26^ < 0.28 for 8q p = 0, but it is still not ruled out because our model 
predicts S^p = —177°. Also, the predicted value r 23 = 0.0352 is larger than the center value 
given in Eq.(5.3) by 3a. Rather, the predicted value is near to the old value r 23 = 0.031 ± 0.005 
given in the second literature in Ref. [8j • We consider that the value of m s is still controversial. 
We can also predict neutrino masses 

m ul ~ 0.0056 eV, m v2 ^ 0.0118 eV, m u3 ~ 0.0507 eV, (5.11) 

by using the input value [12j Am 2 2 ^ 0.00243 eV 2 . We also predict the effective Majorana 
neutrino mass [13] (m) in the neutrinoless double beta decay as 

(m) = {rmU^ + m 2 Ul 2 + m 3 C/ 2 3 | ~ 3.2 x 10" 4 eV. (5.12) 

6 Concluding remarks 

In conclusion, by assuming U(3)xS3 family symmetries, we have proposed a new yukawaon 

1/2 

model in which not only the quark mass matrices M u and Md but also the charged lepton 
mass matrix M e is given by a form (<I>o)(l + a /^)( < ^o) (3*o is a fundamental yukawaon) as shown 
in Eqs.(4.9) - (4.14). The model have only 7 parameters for 15 observables. However, we have 
obtained reasonable predictions. 

Since all the mass matrices for quarks and charged leptons are given by the same form 
($o)(l + a /^0(^o)> we can consider a possibility that if we define (Yo( a /)) = ( < ^ > o)(l + a/^)( < ^o) 
by introducing a yukawaon Yq, all the effective Yukawa coupling constants can be given by 
a single yukawaon VEV (Yq («/))■ However, this idea runs into a stone wall, because we cannot 
distinguish Yo(a,d) from Yo(a e ) in the expression of (Yr), Eq.(4.14). Besides, the VEV of the 
"single" yukawaon Yq is dependent on the parameter o/. This is contradictory to the idea in 
the original yukawaon model that Yj ^ is given by a single yukawaon VEV without including a 
free parameter. Therefore, in the present paper, we have considered that the yukawaons Yj are 
distinguished by the fields 5/ in ^qSj^o as shown in Eq.(2.8). 

One of motivations of the yukawaon model [1] was to understand a charged lepton mass 
relation [J3] m e + m M + m T = (2/3) (^/m^ + ^/m^ + ^/m T ) 2 by considering (Y e ) oc ($> e )(Qe)- 
However, in the present model, since (Y e ) is given by (#o)(l + a eX){&o), the scenario for the 
charged lepton mass relation must be abandoned, and we have to search an alternative scenario 
for the charged lepton mass relation. Nevertheless, apart from this problem, it seems that the 
present new yukawaon model offers us a promising hint for a unified mass matrix model for 
quarks and leptons. 
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